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THE MIXING TIME FOR SIMPLE EXCLUSION 



By Ben Morris 

University of California, Davis 

We obtain a tight bound of 0(L 2 logfc) for the mixing time of 
the exclusion process in Z d /LZ d with k < |L d particles. Previously 
the best bound, based on the log Sobolev constant determined by 
Yau, was not tight for small k. When dependence on the dimension 
d is considered, our bounds are an improvement for all k. We also 
get bounds for the relaxation time that are lower order in d than 
previous estimates: our bound of 0(L 2 logd) improves on the earlier 
bound 0(L 2 d) obtained by Quastel. Our proof is based on an aux- 
iliary Markov chain we call the chameleon process, which may be of 
independent interest. 



1. Introduction. Let G = (V, E) be a finite, connected graph and define a 
configuration as follows. In a configuration, each vertex in V contains either 
a black ball or a white ball (where balls of the same color are indistinguish- 
able), and the number of black balls is at most \V\/2. The exclusion process 
on G is the following continuous-time Markov process on configurations. For 
each edge e at rate 1, switch the balls at the endpoints of e. 

Note that since the exclusion process is irreducible and has symmetric 
transition rates, the uniform distribution is stationary. Let C denote the 
space of configurations and for probability distributions fj,,v on C, let 

\\n-v\\=maxn(Q)-u(Q)= min P(X ^Y) 

QCC X~u,Y^i/ 

be the total variation distance. Denote the mixing time by 

(1) r mix = ml{t : \\K\x, ■) -U\\< \ for all x G C}, 

where K is the transition kernel for the exclusion process and IA is the 
uniform distribution over configurations. 
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This paper is concerned with bounding the mixing time in the important 
special case where G is the d-dimensional torus Z rf /LZ rf ; we call this process 
simple exclusion. 

The exclusion process is a widely studied Markov chain, with connections 
to card shuffling [1, 14], statistical mechanics [2, 6, 11, 13] and a variety 
of other processes (see, e.g., [5, 8]); it has been one of the major examples 
driving the study of mixing times for Markov chains (see, e.g., [1, 3, 5, 14]). 

Quastel [11] and a subsequent independent paper by Diaconis and Saloff- 
Coste [3] bounded the spectral gap for simple exclusion by comparing the 
process with Bernoulli-Laplace diffusion (i.e., the exclusion process on the 
complete graph), obtaining a bound of 0(L 2 d) for the relaxation time (i.e., 
inverse spectral gap); by a standard theorem (see, e.g., [12]), this implies 

Tmix < c[L 2 d] log 

where k is the number of black balls [so that (\ ) is the number of states in 
the chain] and c is a universal constant. 

An improved dependence on the number of states can be obtained using 
the log Sobolev constant. In a famous result of Yau [15], which built on the 
earlier work of Lu and Yau [9], the log Sobolev constant for simple exclusion 
was shown to be 0{cdL 2 ) : where Cd is a constant that depends only on d; 
by a standard theorem (see, e.g., [12]), this implies 

(2) r mix < c d L 2 log log 

(An analysis that combines the comparison technique [3, 4] with Lee and 
Yau's [7] estimate for the log Sobolev constant for Bernoulli-Laplace dif- 
fusion shows that we can take Cd = cd for a universal constant c when 
k = L d /2.) 

The bound (2) was the best known, but it was not tight for small k; in this 
case, the log Sobolev constant, which is based on a single analytic inequality, 
does not capture the faster mixing of the process in its early stages. 

The main contribution of this paper is to give a bound of the correct 
order, namely 0(L 2 log k), for the mixing time. We accomplish this using 
an auxiliary Markov chain called the chameleon process, a variant of the 
exclusion process in which the balls can change color. Using a relationship 
between the transition probabilities of the two processes (see Lemma 1 be- 
low), we give a detailed analysis of the mixing during the various stages and 
obtain a tight bound. 

As an added benefit, when dependence on the dimension d is considered, 
our bounds are an improvement for all values of k, both for the mixing time 
and for the relaxation time. Table 1 shows the improvement when k = L d /2. 
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We have included the special case G = {0, 1} above to illustrate the 
dependence on d. In this case, Wilson [14] showed that the mixing time is 
at least of order d (see [14]), so our bound is within a factor 0(logd) of the 
true answer. Diaconis and Saloff-Coste [3] wondered if the factor d in the 
estimate r re i ax < cd was necessary. We show that it can be replaced by log d. 

We note that the analyses of Quastel, Diaconis and Saloff-Coste and Lee 
and Yau apply to the box {0, 1, ... , L} d in addition to the torus. Our results 
also hold for this and other variants, but in the present paper we will only 
treat the case where G is the torus. 

There is a connection between our key tool — the chameleon process — and 
evolving sets, which were used by the author and Peres [10] to prove heat 
kernel bounds based on isoperimetric inequalities. The reader who is familiar 
with [10] might recognize the chameleon process as a variant of an evolving 
set process (in somewhat disguised form). This connection will not be used 
explicitly. 

2. Interchange and chameleon processes. 

2.1. Interchange process. We will now describe a related Markov chain, 
the interchange process, which will be useful. Let n = \V\. In a configuration 
of the interchange process, balls labeled 1,2, ... , n occupy the vertices of G. 
Each vertex contains exactly one ball, which is colored either black, white, 
red or pink. The transition rule is the following: For each edge e at rate 2, 
flip a coin. If the the coin lands heads, switch the balls at the endpoints of 
e; else do nothing. More precisely, there are independent Poisson processes 
of rate 2 for every edge; when the Poisson process corresponding to edge e 
has an arrival, the balls at the endpoints of e are switched w.p.1/2. When 
there is an arrival that corresponds to e, we will say the clock of edge e rings 
or simply edge e rings. 

For a ball i £ {1, . . . ,n}, define the trajectory of i as the map t — > vt(i), 
where vt(i) is the position of ball i at time t. Note that since each pair of 
neighboring balls is switched at rate 1, the unordered set {vt(l), ■ ■ ■ ,v t (k)} 
behaves like the set of black balls in the exclusion process. 

Table 1 







''"mix 


''"relax 


G 


Old bound 


New bound 


Old bound New bound 


Z d /LZ d 


cd 2 L 2 log(L) 


cd log(d)L 2 log(L) 


cdL 2 clog(d)L 2 


{0,l} d 


cd 2 


crflog(rf) 


cd clog(d) 
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2.2. Chameleon process. The chameleon process is a variant of the in- 
terchange process in which the balls can change color. In the chameleon 
process, the balls move in the same way as in the interchange process; how- 
ever, in addition, each time the clock of an edge e rings a two-step recoloring 
operation is performed. The first step, which we call pinkening, takes place 
when the balls at the endpoints of e are red and white, respectively; in this 
case they are both recolored pink. The second step, which we call depinking, 
takes place when there are a large number of pink balls; in this case all of 
the pink balls are collectively recolored either red or white with probabil- 
ity 1/2 each. More precisely, the transition rule is the following: When the 
clock of edge e rings, switch the balls at the endpoints of e if and only if 
the interchange process switches them. In addition, perform the following 
two-step recoloring operation: 

1. (Pinkening step) If the endpoints of e contain a red and a white ball, 
respectively, then recolor both balls pink. 

2. (Depinking step) If the number of pink balls (including the pink balls 
created in step 1) is at least as large as either 

• the number of red balls or 

• the number of white balls, 

then flip a coin. If the coin lands heads, recolor all pink balls red; else 
recolor all pink balls white. 

Starting colors. We will always start with the balls colored in the fol- 
lowing way. For some b S {0, . . . , k — 1}, balls 1, . . . , b are colored black, ball 
b + 1 is colored red and the remaining balls are colored white. Note that the 
number of black balls b remains fixed for the duration of the process, and 
that the process is symmetric with respect to the roles of the red and white 
balls. 

2.3. Overview of the chameleon process. We now introduce some nota- 
tion. For a configuration M of the chameleon process, let R(M), W(M) and 
P(M) denote the set of balls in M that are red, white and pink, respectively, 
and let 

r{M) = \R{M)\; w{M) = \W{M)\; p(M) = \P(M)\. 

For a chameleon process {Mt,t > 0}, define Rt := R(Mt) and so on. 

We will call the times at which the chameleon process recolors its pink 
balls depinking times. Let Tq = and for j £ {1, 2, . . . }, let Tj be the jth 
depinking time. We will now describe the behavior of the chameleon process 
between time Tj and TV+i- At time Tj, there are only black, white and red 
balls (since depinking changes all pink balls to either red or white). Let b be 
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the number of black balls and let m = n — b. Assume that there are at least as 
many white balls as red balls at time Tj , so that < ^ (the behavior can 
be described symmetrically when > y). From time Tj until time 2y+i, 
the process behaves like the interchange process, except that when an edge 
connecting a red ball to a white ball rings, the two balls are each recolored 
pink (in the pinkening step). Pink balls accumulate until there are at least 
as many pink balls as red balls, and then the next depinking occurs at time 
Tj + \ (at which point the pink balls are recolored either red or white). Since 
each pinkening between time Tj and time Tj + \ increases the number of pink 
balls by two and decreases the number of red balls by one, the depinking 
occurs when the number of red balls falls below two-thirds of its value at 
time Tj. Hence we have 

r Tj+1 = r Tj ± 

When there are more red balls than white balls at time Tj, the process 
behaves as described above, with the roles of red and white balls reversed. 
Hence, in general we have 

TT j+1 =rr j ±A(r T] ), 

where A(x) := |" a:A ( ? ^~ a: ) "| _ The decision whether to choose + or — is made 
by a coin flip that is independent of the trajectories of the balls. Hence, the 
process {r^ - J ' > 0} is a Markov chain. (Note that it is also a martingale.) 

3. Main theorem. Define the redness p(u, M) of vertex u in configuration 
M by 

!1, if the ball at u is red; 

1/2, if the ball at u is pink; 
0, otherwise. 

If {M t ,t > 0} is a chameleon process, define pt(u) := p(u,M t ). The quantity 
X) u Pt{u) =r t + \pt is a measure of the amount of "red paint" in the system 
at time t, if we think of the pink balls as being painted with a mixture of 
red and white paint. Similarly, the quantity wt + \pt measures the amount 
of white paint. The Markov chain : j = 0, 1, . . . } has absorbing states at 
and m. Let A be the event that it absorbs at m, the system eventually 
consisting of only red and black balls with all "white paint" gone. Since {rj^ } 
is a martingale and rj- = tq = 1, the optional stopping theorem implies 
P(A) = Note that A is completely determined by the coin flips that 
are performed during the depinking steps and hence is independent of the 
trajectories of the balls. 

Define v t (l, . . . ,k) := (v t (l), . . . , v t (k)) and let X t = Vf(l, ...,&) be the vec- 
tor of positions of balls 1, . . . , k at time t. Since the unordered set {vt(l), • • • , 
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vt(k)} behaves like the set of black balls in the exclusion process, any up- 
per bound on the mixing time of the process {X t , t > 0} also applies to the 
/c-particle exclusion process. 

Let nt denote the distribution of Xt. Then p^^U as t — > oo, where U is 
the distribution of k uniform samples without replacement from V. 

For j <k, we write 

//t(ui,...,Uj) :=P(v t (l,...,j) = (ux,...,Uj)); 

(3) 

Ht(v,j\ui, . . . ,Uj-i) :=P(v t (j) =Uj\v t (l,...,j - 1) = (ui,...,Uj-i)), 

with analogous notation for IA. Note that each of the conditional distribu- 
tions Ht{-\ui, . . . ,Uj—i) converges to uniform as t — > oo. When all of these 
distributions are "very close" to uniform, then pt is close to IA. (In fact, we 
only need the conditional distributions to be close, "on average"; Lemma 12 
in Section 5.3 formalizes this.) Lemma 2 below relates the mixing of each of 
the conditional distributions pt(-\ui, . . . ,u b ) to a corresponding chameleon 
process with b black balls. It gives a bound based on the expected amount 
of "white paint" in the system at time t, conditional on A. 

To prove Lemma 2, we will need the following lemma, which indicates the 
fundamental relationship between the interchange process and the chameleon 
process. Recall that pt(u) denotes the redness of vertex u at time t. 

Lemma 1. Fix b G {0, ... ,k — 1} and let zi, . . . , 2^ + i be a sequence of 
distinct vertices in V . Consider a chameleon process with b black balls. Then 

(4) pt(z b+ i\zi, ...,z b ) = E(p t (z b+1 )\v t (l, ...,&) = (zi,.. -,Zh)). 

Proof. Let n < ti < ■ • ■ be the times when edges ring. For j > 1, let 
{Ml , t > 0} be the process that behaves like the chameleon process before 
time Tj and behaves like the interchange process at time t 3 - and afterward 
(i.e., all recoloring is suppressed from time Tj onward). Note that {M/ : t > 0} 
is the interchange process. We claim that for every j > 1 we have 

(5) p t (z b+ i\zi,...,z b ) = E(p(z b+ i,Ml)\v t (l, . . . ,b) = (zi,...,z b )), 

where p(z b +i, Mf) denotes the redness of vertex z b+ i in configuration M/. 
This implies (4) because Mt is the a.s. limit of M( as j — > oo. We will show 
that (5) holds for all j > 1 by induction. The base case j = 1 is trivial because 
{M/} is the interchange process and has no recoloring, so ball b + 1 is always 
the only red ball (and there are never any pink balls). Thus, it remains to 
show that for every j > 1 we have 



(6) 



E(p(z b+ i,Mi)\D) = E(p(z b+ i,Mt l )\D), 
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where D is the event that vt(l, ■ ■ ■ , b) = (z\ , . . . , z b ) . To show this, we intro- 
duce a new process {Mj + } that behaves like {Ml } except that it performs 
a pinkening step at time Tj+\ (but no depinking then). We first show that 

(7) B(p(z b+1 ,Mi)\D) = B(p(z b+1 ,Mi + )\D). 

Let ej+i be the edge that rings at time Tj + % and let {Ml } be the process 
that behaves exactly like {Ml}, except that if the balls at the endpoints of 
ej+i are red and white, it switches them if and only if {Ml } does not switch 
them. Clearly Ml has the same distribution as Ml and hence 

(8) E(p(z b+1 ,Mi)\D)=E(p(z b+1 ,Mi)\D). 
But note that 

\p{z b+1 ,Mi) + lp(z b+1 ,Mi) =p{z b+1 ,Mi + ). 

Taking a conditional expectation given D and combining with (8) gives (7). 
It remains to show that 

(9) E(p(z b+1 ,Mi + )\D) = B(p(z b+1 ,Mi +1 )\D). 

Recall that M J+1 behaves exactly like M J+ except that it performs a de- 
pinking step at time Tj + \. Let M^ +l be the process that behaves exactly like 
AP +1 , except that if it depinks at time Tj+i then it makes the opposite recol- 
oring choice (i.e., if M^ +1 recolors its pink balls white, then M J+1 recolors 
its pink balls red and vice versa). Clearly M^ +1 has the same distribution 
as M J+1 and hence 

(10) E(p(z b+1 ,Mi +1 )\D) = E(p(z b+1 ,Mi +1 )\D). 
But note that 

\p{z b+1 M +1 ) + lp(z b+1 ,Mi +1 )=p(z b+1 ,Mi + ). 

Taking a conditional expectation given D and combining with (10) gives 
(9). □ 

We are now ready for Lemma 2. Let U(-\z\, . . . , z b ) denote the uniform 
distribution over V — {z±, ... ,z b ). 

Lemma 2. Fix t > and let (Z±, . . . , Z b ) = vt(l, . . . ,b). Consider a cha- 
meleon process with b black balls. Then 

(11) E(\\p t (-\Z l ,...,Z b )-U(-\Z 1 ,...,Z b )\\)<-E(w t + l P X 
where \\ ■ \\ denotes variation distance and we write E(-) for E(-|^4). 
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Note that since Zi,...,Z b are random variables (whose joint distribution 
is governed by p t ), so is the quantity ||^(-|Zi, . . . , Z b ) —U(-\Zi, . . . , Z b )\\ ; the 
left-hand side of (11) is the expectation of this random variable. 

Proof of Lemma 2. Recall that for probability distributions v\,V2, 
the total variation distance 

ll^i - v 2 \\ = sup(i^ 2 (Q) - v x {Q)). 
Q 

Fix distinct vertices z±, . . . , z b and a subset Q C V — {z±, . . . , z b }. Let 

z = (zi,...,z b ); Z = (Z 1 ,...,Z b ) 
and m = n — b. Lemma 1 implies that 

(12) U(Q\ Zl , ...,z b )- m(Q\ Zl , . . . ,z b ) = M - E( Pt (Q)\Z = z), 

m 

where we write pt(Q) for J2xeQ Pt( x )- But 

B(p t (Q)\Z = z)>-E(p t (Q)\Z = z), 
m 

since P(.A) = — ■ Hence the right-hand side of (12) is at most 

-E(\Q\-p t (Q)\Z = z) = -E(\QnWt\ + hQnP t \\Z = z 

< ^E\w t + ^p t \Z = z 

Since this holds for all Q, we have 

\\p, t (-\zi,...,z b ) -U(-\zi,...,z b )\\ < —lt(wt + -pt\Z = z) 

m \ 2 J 

and hence 

(13) E(||/x t (-|Z 1 ,... ) Z 6 )-W(-|Zi,... ) Z 6 )||)<^E('EL t + ^z)Y 

Recall that the event A is independent of the trajectories of the balls (and 
hence independent of Z). It follows that we can rewrite the RHS of (13) as 

For a chameleon configuration M, let S(M) = (r(M) + ^p(M)) and let 

s(M) = For a configuration M, let W(M) = E A/ (Ti), where T x is 

the first depinking time. So W< is the expected waiting time until the next 
depinking after time t, given Mt- In Section 4, we show that 

m<g{S t ), 



7 = 

Note that < 7 < y/2 - 1. Let 
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where g(S) = clogd(S + b) 2 l d for a universal constant c. Let 

_ 2 =-iY 

4/3 + v/2/3 / 

[ ' 9(S(M)) 
Then 1 < Y < 1 + 7. For < u < 1 define it" = min(it, 1 — u). We have 

E(l - a t ) < E(l - at )Y t 



where the second inequality can be seen by considering the cases st > \ and 
s t <\ separately. Define Z(M) = ^/ffV (M) and let L t = E(Z t ), so that 
E(l - s t ) < L t . Note that 

L t = J2 P i M t = M)Z(M) 

M 

is differ entiable. Define / by 



\g{{l +1 ) 2 m / z 2 y\ ifz>V2, 
i(clogd) _1 L~ 2 , otherwise. 



Note that g is increasing and, for z > \/2, we have g((l + i y) 2 m/z 2 ) < g(m) = 
clogdL 2 [recall that (7 + 1) 2 < 2]. It follows that / is increasing. For Z t < y/2, 
we have 

f(Z t ) = Kclogd)- 1 ^ 2 < l(clogd(S t + b)y 2,d = \g(S t )-\ 
and, when Zt > \/2, we have st = s\ and hence 

m)=\mfs t )- l <\ g{St r>. 

Thus /(Z t ) < for all t. 

Lemma 3. The differential equation 
(14) L[ < - 7 L t f(L t /2) 

is satisfied. 
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Proof. First, we will show that for every t > 0, 

(15) E(Z t+£ \M t ) < Z t (l - 2rfef(Z t )) + 0(e 2 ). 

By the Markov property, it is enough to show that for any starting config- 
uration Mq we have 

(16) E Mo (Z e ) < Z (l - 2 7 e/(Z )) + 0(e 2 ). 

Note that St is constant on each interval [Tj, and that st^ = for 

all j. It follows that {st :i > 0} is a martingale and hence P(A|so = s ) = — 
for every s. Using this, it is easy to check that 

(17) E(h(M t ))=E^h(M t )^j 
for every function h and t > 0. Thus 

(18) E(Z £ ) = E(Z £ ^) =lE( V / 4y e ). 
We claim that 

(19) V(sfiY £ ) <(l- ^ZL^ y/j Q Y + 0(e 2 ). 
Combining this with (18) gives 



<(l-2 1 ef(Z ))Z + O(e 2 ), 

which verifies (16). 

It remains to verify (19). Let {M t :t > 0} be the process that behaves like 
Mt, except that it does not perform any depinking steps. So Mt pinkens, but 
does not depink and hence accumulates pink balls. Let p(-), f (•) and w(-) out- 
put the number of pink, red and white balls, respectively, in a configuration 
M, and define 

VW(M) =E A7 (inf {t:pt >mm(f t ,w t )}), 

where we write E^(-) := E(-|M = M). Let T x be the first time that M t 
depinks. If we start with Mq = Mq, then for every t > we have 



W(M, 



(W(M t ), if t<T i; 
\0, ift>Ti. 
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By analogy with the previous discussion, let S = (w + \p), let s = and 
define Y(M) 

min(f, w). We will show that 



1 + ■ Fix a starting configuration Mq = Mq with p < 



(20) 

and 
(21) 



E(Vfc) < 1 



2g(S ) 



lsl% + 0(e 2 ) 



E( v / 57y £ )<E(v^y £ ) + o( £ 2 



These will imply (19), since Yq = Yq and so — sq. Note that since Mt never 
depinks, st = s is constant and we have 



(22) 



E(Vsly e )-E(V4^c 



^ 7 E(W(M e )-W(M )) 



Let B = {M :p > min(f, w)}. For any configuration M, the quantity W(M) 
is the expected waiting time to hit B starting from M. Hence Proposition 
11 below implies that E(W(M e ) - W(M )) = -e + 0(e 2 ). Combining this 
with the fact that Yq < 2 gives 



(23) 



E( V / s|y e )-E(V4 i o)< vs 



\Y + O{e 2 ) 



MS) 



and (20) is verified. It remains to verify (21). Let D be the event that there 
is exactly one transition in [0, e] and let F be the event that M £ G B. Note 
that P(FuD)>l-0(e 2 ), and on the event F c we have M £ = M £ and 
Y e = Y e . Thus, it is enough to show that 



(24) 



E(\f slY £ \D <1 F) <E(\/ slY e \D D F) 



On the event D n F, M £ is obtained from M £ by performing a depinking 
operation. Fix M £ B with s < ^ and let M be the configuration obtained 
from M by depinking. Then M is a random configuration with two possible 
values, since the pink balls can be recolored either red or white. Let Y = 
Y(M) and s = s(M). Since Y < 1 + 7 and s" < s, we have 

E(v / sly) < (l+ 7 )E(y / i) 



= (l+7)(^ V /F+A(5) + i v /s 
<(l+7)(- 



A(S) 



/I73 + v/2/3 



where the second inequality holds because A(s) > s/3. But, s = and 
y > 1, so E(v / s«y) < vW. This holds for all M e B with s < 1/2. By 
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reversing the roles of red and white balls, we can argue similarly in the case 
s > i. This verifies (24), which in turn verifies (15). 

Taking the expectation (conditional on A) of both sides of (15), we get 

(25) L t+£ <L t - 2 7 eB(Z t f(Z t )) + 0(e 2 ). 
We will use the following lemma from [10]. 

Lemma 4 ([10]). Suppose that Z>0 is a nonnegative random variable 
and f is a nonnegative increasing function. Then ~E(Zf(2Z)) > -s- • f(EZ). 

Lemma 4 implies that B(Z t f(Z t )) > ^ff(^f)- Plugging this into (25), we 
get 

(26) L t+£ <L t - e 7 L t f(L t /2) + 0(e 2 ) 
and the proof is complete. □ 

Remark. Note that since / is bounded below by ^((logtiL 2 )" 1 ), Lem- 
mas 2, 3 and 12 imply that the spectral gap for the simple exclusion process 
must also be at least of order (log dL 2 )" 1 . 

We are now ready to prove our final ingredient, namely, a bound on the 
RHS of (11) when G is the torus: 

Lemma 5. Let G = Z d /LZ d . Fix b € {0, . . . , k — 1} and consider a chameleon 
process on G with b black balls. There is universal constant C such that if 
t>C[d\ogdL 2 + \ogdL 2 \ogk], then 

(27) Ig^ + lp,)^. 

Proof. We will need the following lemma from [10]. 

Lemma 6 ([10]). Suppose that Lt is a nonnegative, differ entiable func- 
tion oft that satisfies the differential equation L' t < —jLtf(Lt/2). Then for 
every 

rio/2 fa 

t> 

we have L+ < 5. 



s/2 -yzf(z)' 



Recall that the number of vertices n is L d . Since L = Z = ^=^o < 2^pa, 
Lemmas 6 and 3 imply that Lt < -h- for 

" dz 



(28) t > 



i/sk jzf(z) 
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We need to show that this integral is at most 
(29) Cd log dL 2 + C log dL 2 log k 

for a constant C. We can write the integral as 



(30) 



dz 



+ 



dz 



'l/8fc 7^/(^) 7*/(z) 

. boun 
when z < y/2. 



+ 



dc 



i/6 jzf(z) 



h + h + h- 



We will bound the terms Ix,l2 an d I3 separately. Since = 4c7 1 log dL 2 



I x = 4 C1 ~ 1 bg dL 2 log(8v / 2it) < a log dL 2 log A; 

for a constant a and hence the second term in (29) dominates I\ if C is 
sufficiently large. To bound I2 and I3, note that for z > \/2 we have 



(31) 



1 



■ 4c7 1 log d 



for a constant (3. Hence 



m(l+7) 2 



2/d 



+ b\ < (3 log d • max [ —r , 6 



2/d 



'm/b 



(3 log d • - dz 



2/dj 



"( 1+4 / d ) dz). 



< 0\ogdm 2 / d ( 

But, m 2 / 6 ' < n 2 / d = I? and an elementary calculation shows that the integral 
in the second line is at most d. Hence the first term in (29) dominates I2 if 
C is sufficiently large. It remains to bound I3. Using (31) we get 



(32) 



But since 



/v n 1 
_f3\ogd-b 2 ' d -dz 
z 

< (3log db 2/d 



-log(n/m) + ^logfe 



m > n/2; 6 2/d < n 2/d = L 2 ; 6 < fc, 



we have I3 = 0(1? logdlog k). Thus I3 is dominated by the second term in 
(29) if C is sufficiently large. 

We have shown that I\, I2 and I3 are all of the correct form. This com- 
pletes the proof of Lemma 5. □ 



We are now ready for our main theorem. 
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Theorem 7. The mixing time r m j x for k-particle exclusion in Z d /LZ d 
satisfies 

(33) r mix < C [d log dL 2 + log dL 2 log k] 
for a universal constant C . 

Proof. Let G = Z d /LZ d . Fix t > and suppose that (Zt, . . . ,Z k ) ~ n t . 
Lemma 12 in Section 5.3 gives 
k-i 

(34) <J2V(\\M-\Zi, . . . ,Z b ) -Ui-lZu . . . ,Z b )\\), 

6=0 

but by Lemmas 2 and 5, each term in the RHS of (34) is at most ^ whenever 
t > C[dlogdL 2 + log dL 2 log /c] for a universal constant C. Hence the sum is 
at most |. □ 

4. Waiting to depink. Let G = Z d /LZ d . Consider a chameleon process 
on G with b black balls. Let T\ be the time of the first depinking. The 
principal result of this section is an upper bound on E(Ti), valid uniformly 
over starting states M . 

Recall that Wt and Rt denote the set of balls that are white and red, 
respectively, at time t. Say an edge is conflicting at time t if one of its 
endpoints contains a ball from Wq and the other contains a ball from Rq. 

We will need the following lemma. 

Lemma 8. There exist universal positive constants c and B, with B < 
2/3, such that if r ~ uniform[0, c\ogd(b + S) 2 / d ], then the expected number 
of edges conflicting at time r is at least Br^d. 

Proof. Let G= (V,E) be a connected graph. Define the random walk 
on G with rate A as the continuous-time Markov chain on V with the fol- 
lowing transition rule: If the current state is x S V, traverse each edge e 
incident to x at rate A. Let G = Z d /LZ d and let {Zt :t > 0} be the random 
walk on G with rate 2. Let G' be the subgraph of G induced by removing 
the origin and let {Z' t :t > 0} be the random walk on G' with rate 2. 

Consider the chameleon (or interchange) process on G. Clearly, for any 
distinct pair of balls x and y, we have 

(35) P(v t (x) is adjacent to v t (y)) = P(Z' t £ B\Z' = v (x) - v (y)), 

where B is the set of neighbors of the origin in G. 

The Markov chain {Zt} is easier to work with than {Z^}. They are related 
in the following way: Fix z £ V — {0} and a > 0, and let t ~ uniform [0, a]. 
Then 

(36) P{Z' T eB\Z^ = z)>-p(Z T eB\Zo = z). 
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The reason that (36) holds is that the "wasted time" Z± spends in lowers 
the probability that Z T £ B. To give a formal justification for (36), we will 
show that \{t < a:Z[ £ B}\ stochastically dominates \{t < a: Z t £ B}\. Let 
Z t be the process obtained from Z t by "skipping past" all visits to the origin, 
that is, the process obtained by replacing transitions to by transitions to 
a uniformly chosen element of B. 

Note that in the natural coupling (where the trajectory of {Zt} is con- 
structed from that of {Zt} by simply deleting every visit to the origin), we 
have 

\{t<a:Z t eB}\ > \{t<a:Z t €B}\. 

Also, clearly we can couple (Z t , Z' t ) so that 

{t:Z t eB} = {t-.Z' t eB}. 

It follows that \{t < a: Z[ £ B}\ stochastically dominates \{t < a:Zt £ B}\ 
and (36) is verified. 

Let p(-, •) be the transition kernel for {Zt} and, for e > define 

r(e) = mi{t:p\x,y) < \e Vx,y}. 

We will need the following easy proposition, which is proved in the next 
section. 

Proposition 9. For any e > \ , we have 

T(e)<e~ 2/d Dlogd, 
where D is a universal constant. 

We are now ready to complete our proof of Lemma 8. Let c = 20D, where 
D is the constant that appears in Proposition 9, and let B = 1/72. For 
x £ i?o, let N x = {y £ Wq :v T (y) is adjacent to v T (x)}. It is enough to show 
that for every x £ Rq, we have E(N X ) > Bd. But 

(37) E(N X ) = f ( v t(x) is adjacent to v T (y)) 

y ew 

(38) = Y, P ( Z 'r e B\Z' = v (x) - vo(y)) 

yeW 

(39) > J2 p ( z r e B\Zo = v (x) - v (y)) 

y ew 

(40) =E(p T (A,B)), 
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where A = {vo(x) — Vo(y) :y £ Wo}. The second line follows from (35) and 
the inequality follows from (36), but the symmetry of {Zt} implies that the 
quantity (40) is 

E(p T (B,A)) = \B\P UB (Z T eA) 
> dmmP x (Z T eA), 

xGV 

where Ub denotes the uniform distribution over B. 

We will now complete the proof by showing that ~P X (Z T £ A) > 1/72 for 
all x £ V . Fix x E V . Since r ~ uniform[0,clog(i(6 + S) 2 / d ] and c = 20D, we 
have 



P x (Z T eA) > \V X {Z T eA\r> 10Dlogd(b + S) 



Thus, it is enough to show that p\x,A) > 1/36 for t > WD log d{b + S) 2/d . 
Let A c denote the complement of A in V . Then 

p\x,A) = l-p\x,A c ) 
(41) > 1 - \A C \ ma,xp\x,y) 

y&V 

= 1 — (b + ro +po)maxp*(x,y). 



5n/6. Plugging e = | /, , 1 , — r into Proposition 9 shows that for t > Dlogd- 



Since b < n/2 and w o > r$ > po, we have u>o > n/6 and hence (b + ro +Po) < 

5 1 

6 (b+r +p ) 

[|(& + r- +p )] 2/d , we have 

35 1 



P (z,y) < 



36 (b + ro+po) 



for any y and hence p t (x,A) > 1/36 by (41). Finally, note that S > (ro + 
Po)/2 and hence 

D log d • [f (6 + r + p )] 2/d < -D log d[3(6 + S)] 2 / d 

< lODlogc^ + S) 2 ^. □ 

We are now ready for the main result of this section. Recall that T\ 
denotes the time of the first depinking. 

Lemma 10. Let S = r$ + ^po and suppose S £ {0,m}. Then 
(42) E(Ti) < clogd(6 + S) 2 ^ d 

for a universal constant c. 



THE MIXING TIME FOR SIMPLE EXCLUSION 



17 



Proof. Since the chameleon process is symmetric with respect to the 
roles of the red and white balls, and since the expression on the RHS of (42) 
is increasing in S, we may assume that S < to/2. Note that rt + \pt = S for 
all t G [0,Ti). For (3 > 0, let 

Ta = min{£ : at least ftr^ balls have been pinkened by time t}. 

We will show that there exist universal constants c, (3 > such that 

(43) P(Tp<clogd(b + S) 2 / d )>f3. 

Note that the lemma follows from (43) since we can incorporate an extra 
factor into c. 

It remains to verify (43). Recall that Wt and Rt denote the set of balls that 
are white and red, respectively, at time t. Say edges e and / are neighbors 
if they are incident to a common vertex. Note that for any edge e, 

P(e rings for some time t G [0, 1] and before any of its neighbors) > — . 

5a 

(44) 

To see this, fix an edge e. Let A\ be the event that e rings before any of 
its neighbors and let Ai be the event that e rings for some t G [0, 1] . Since 
no edge in 7i d /L7i d has more than 4d — 2 neighbors, we have P(^4i) > zxrp 
and since the edges ring at rate 2, we have P(^2) = 1 — e ~ 2 ■ Finally, note 
that A\ and A<i are positively correlated and hence 

Recall that an edge is conflicting at time t if one of its endpoints contains 
a ball from Wo and the other contains a ball from Rq. Let N denote the 
number of edges in G that: 

1. Conflict at time r. 

2. Ring before any of their neighbors after time r. 

3. Ring for some t G [r, r + 1]. 

Lemma 8 says that the expected number of edges that conflict at time r 
is at least BrQd for a constant B. Combining this with (44) and using the 
strong Markov property, we get 

Let J be the number of balls pinkened before time 2\ A (r + 1). Note that 
J > N. Thus E(J) > E(iV) > Bro/5. Hence, an application of Markov's in- 
equality to the nonnegative random variable 2ro — J gives 



10 J ~ 20 - B 
B 

> 20" 
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Finally, note that on the event that J > Btq/IO, we have Tb/\q < t + 1. It 
follows that 

P(r B/10 < clog (2(5 + b) 2/d + 1) > B/20. 

Incorporating the (small) additive constant 1/ log (2(5 + b) 2 l d into c and 
setting (3 = B/20 verifies (43) . □ 

5. Technical details. We now prove three technical results that were 
needed in previous sections. 

5.1. Easy proposition used in the proof of Lemma 3. 

Proposition 11. Let Xt be a continuous-time Markov chain on a finite 
state space V and let B CV be a subset of states. For x £ V , let W(x) = 
E x (inf{2 :Xt £ B}) be the expected time to hit B starting from x. Then, for 
every x £V — B , we have 

B X (W(X £ ) - W(x)) =-e + 0{e 2 ). 

Proof. Fix x & V — B and let Xq = x. Let q x be the rate of transitions 
out of x, let Y be the next state visited after x and let W = inf{t : Xt = Y}. 
Then W has an exponential distribution with parameter q x . It follows that 
W(x) = l/q x + E(W(y)). Let D be the event that the chain makes exactly 
one transition in [0, e]. Then 

B(W(X £ ) - W(x)\D) = E(W00) - W(x) = -l/q x 

and P(L>) = eq x + 0(e 2 ). Hence 

B(W(X £ ) - W(x)) = (eq x ) ■ (-l/q x ) + 0(e 2 ) 

= -e + 0{e 2 ). □ 

5.2. Probability bounds for random walk on the torus. 

Proof of Proposition 9. We adopt the notation of the previous sec- 
tion. Let Xt be the random walk of rate 2 on the cycle Cx of length L. It is 
well known that for all j < L, I > and t> j 2 l, we have 

P(^=y)<-(i + e"^) 

J 

for every y £ Cl, where 7 > is a universal constant. Thus the proposition 
holds when (2 = 1 (provided that we replace log d by 1 in the bound). 
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When d > 1, the coordinates of Zt are independent copies of Xf. Thus, 
for every z £ V and t> j 2 l, we have 

V{Z t = z)<j d {l + e^ l ) d . 
Plugging in I = ^ log <i we get 

1 / 1 \ d 1 7 

where we have used the inequality (1 + 1/ d A ) d < |, valid for d>2. The above 
shows that for all j < n, we have 

T (?) £i 1 logd 

Setting e = and D = 4/7 we get 

T(e)<e- 2/d Dlogd. □ 

5.3. Lemma about total variation distance. 

Lemma 12. Fix a graph G, let [i and v be arbitrary probability dis- 
tributions supported on k-tuples of distinct vertices of G and suppose that 
(Zi,.. .,Z k ) ~ (i. Then 

k-l 

(45) \\ f i-u\\<Y,H\H-\Zi,...,Z l )-u(-\Z 1 ,...,Z l )\\). 

1=0 

Proof. Recall that for probability distributions p, and V, the total vari- 
ation distance 

(46) \\j2-V\\= min ^P(W 1 ^W 2 ). 

Thus, for every I and z±,...,Zi, we can construct W\ ~ /j,(-\zi, . . . ,zi) and 
W2 ~ v(-\zi, . . . , z{) such that 

P(W 1 ^W 2 ) = \\ f i(-\z 1 ,...,z l )-v(-\z 1 ,...,z l )\\. 

We couple Z ~ fi with y ~ f as follows. Choose (Yi,Z\) according to the op- 
timal coupling [i.e., a coupling that achieves the minimum in the RHS of (46)] 
and subsequently, for all I with 1 < I < k — 1: if (Z±, . . . , Z{) = (Yi, ...,!/), 
then choose (Zi + \, Y1+1) according to the optimal coupling of fi(-\Z\, . . . , Zj) 
and v(-\Zi, . . . , Z/); else couple (YJ+i, Zi + \) in an arbitrary way. 
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Note that 

fc-i 

(47) P(Z ±Y) = Y d P((^l, ...,Z l ) = (Y 1 ,...,Y l ),Z l+1 ^ Y l+1 ), 

1=0 

but on the event that (Zi, . . . , Z{) = (Y\, ...,Y{), the pair (Zi+±, ij+i) is cho- 
sen according to the optimal coupling of fx(-\Zi, . . . , Zi) and u(-\Z\, . . . , Z\). 
Hence the RHS of (47) is at most 

fc-i 

^Edi/iOiZ!,...,^)-^-^!.---.^)!!)- n 

1=0 u 

6. Lower bounds. In this section, we prove a matching lower bound 
of 0(L 2 logA;) for the mixing time of the exclusion process in Z d /LZ d with 
k black balls. We will accomplish this by showing that there is a universal 
constant c and an exchangeable event A, such that 

U{A)-m{A)>\ 

for every t < cl? log/c, where U and \xt are defined as in Section 3. 

The idea behind the proof is the following: If the black balls start out with 
their first coordinate close to L/4, then with high probability most will still 
be in the left half of the torus at time t = Q(L 2 logk). For convenience, we 
assume that L is a multiple of 8. The proof for the general case is similar. 
Suppose that balls 1, . . . , k are black and the remainder are white. Suppose 
that vt(x) = (v}(x), . . . ,vf(x)), and let the starting configuration be chosen 
to minimize the quantity 

k 

(48) £ 

x=l 

and such that Vq(x) < L/2 for every black ball x (so that the first coor- 
dinate takes the value in preference to L/2). For v £ V, suppose that 
v = (v 1 ,. . . , v d ) and let 

S=L eV :^<v 1 <^y, U=( [ vGV:^<v 1 <Ly 

Note that for every ball x, the process {v\ (x) : t > 0} is just the continuous- 
time random walk of rate 2 on the cycle Cl of length L. Hence P(vt(x) £ 
U)< \ for every black ball x and t > 0. It is also well known that if vo(x) £ 
S, then we have P(vt(x) E U) < ~ — <5exp(— jt/L 2 ) for universal constants 
5,7 > 0. Let 

N t = \{l<x<k:v t (x)£U}\ 



It \ ^ 
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be the number of black balls in U at time t. Since at least half of the black 
balls must start in S, we have E(Aj) < |[1 — <5exp(— ^t/L 2 )] for every t. 

For vertices v £ U , let X t {v) be the indicator for the event that v is 
occupied by a black ball at time t, so that N t = ^ v( zu X l {v). Liggett [8] 
proved that for u^v, the random variables X t (u) and X l {v) are negatively 
correlated, that is, Y l {X t {u)X t {v)) < E(X*(u))E(A*(t>)). It follows that 

E(iVf) = E ( £ X\v) + £ X* („) ] 

<e(£x»] +E E ( x *W) E ( xt W) 

< E(iV t ) + (E(A^j)) 2 

and hence var(A r ^) < E(A^) < k/2. Thus Chebyshev's inequality gives 

nf.r k \ vav(N t ) 
P [Nt > - < 



2 / " ((^/2)e-^ 2 )2 

1 _ 2e 2 ^/ L2 

(M 2 /2)e- 2 Wi 2 ~ 5 2 A: ' 



< 



It follows that for t < L 2 logA;/47 we have P(iV t >k/2)< 2/S 2 Vk, which is 
less than | if k > 64/<5 4 . But, a uniformly chosen configuration will have at 
least k/2 black balls in U with probability at least ^. This completes the 
proof. 

Acknowledgments. We are grateful to A. Bandyopadhyay, N. Berger, 
M. Liu, R. Lyons, F. Martinelli, E. Mossel, Y. Peres, A. Sinclair and P. Tetali 
for useful discussions. 

REFERENCES 

[1] Benjamini, I., Berger, N., Hoffman, C. and Mossel, E. (2005). Mixing times of 

the biased card shuffling and the asymmetric exclusion process. Trans. Amer. 

Math. Soc. 357 3013-3029. MR2135733 
[2] Cancrini, N. and Martinelli, F. (2000). On the spectral gap of Kawasaki dynamics 

under a mixing condition revisited. J. Math. Phys. 41 1391-1423. MR1757965 
[3] DlACONlS, P. and Saloff-Coste, L. (1993). Comparison theorems for reversible 

Markov chains. Ann. Appl. Probab. 3 696-730. MR1233621 
[4] DlACONlS, P. and Saloff-Coste, L. (1996). Logarithmic Sobolev inequalities for 

finite Markov chains. Ann. Appl. Probab. 6 695-750. MR1410112 
[5] Fill, J. (1991). Eigenvalue bounds on convergence to stationarity for nonreversible 

Markov chains with an application to the exclusion processes. Ann. Appl. Probab. 

1 62-87. MR1097464 



22 



B. MORRIS 



[6] Kipnis, C, Olla, S. and Varadhan, S. (1989). Hydrodynamics and large deviations 
for simple exclusion processes. Comm. Pure Appl. Math. 42 115-137. MR978701 

[7] Lee, T.-Y. and Yau, H.-T. (1998). Logarithmic Sobolev inequality for some models 
of random walk. Ann. Probab. 26 1855-1873. MR1675008 

[8] Liggett, T. M. (1985). Interacting Particle Systems. Springer, New York. 
MR776231 

[9] Lu, S. and Yau, H.-T. (1993). Spectral gap and logarithmic Sobolev inequality for 

Kawasaki and Glauber dynamics. Comm. Math. Phys. 156 399-433. MR1233852 
[10] MORRIS, B. and Peres, Y. (2006). Evolving sets, mixing and heat kernel bounds. 

Probab. Theory Related Fields. To appear. 
[11] Quastel, J. (1992). Diffusion of color in the simple exclusion process. Comm. Pure 

Appl. Math. 45 623-679. MR1162368 
[12] Saloff-Coste, L. (1997). Lectures on finite Markov chains. Lectures on Probability 

Theory and Statistics. Lecture Notes in Math. 1665 301-413. Springer, Berlin. 

MR1490046 

[13] Thomas, L. E. (1980). Quantum Heisenberg ferromagnets and stochastic exclusion 

processes. J. Math. Phys. 21 1921-1924. MR575630 
[14] Wilson, D. (2004). Mixing times of lozenge tiling and card shuffling Markov chains. 

Ann. Appl. Probab. 14 274-325. MR2023023 
[15] Yau, H.-T. (1997). Logarithmic Sobolev inequality for generalized simple exclusion 

processes. Probab. Theory Related Fields 109 507-538. MR1483598 

Mathematics Department 
University of California, Davis 
1 Shields Avenue 
Davis, California 95616 
USA 

E-MAIL: morris@math.ucdavis.edu 



